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ABSTRACT 

As part of our ongoing effort to investigate transit timing variations (TTVs) of known 
exoplanets, we monitored transits of the four exoplanets HAT-P-18b, HAT-P-19b, 
HAT-P-27b/WASP-40b and WASP-21b. All of them are suspected to show TTVs 
due to the known properties of their host systems based on the respective discovery 
papers. During the past three years 42 transit observations were carried out, mostly 
using telescopes of the Young Exoplanet Transit Initiative. The analyses are used to 
refine the systems orbital parameters. In all cases we found no hints for significant 
TTVs, or changes in the system parameters inclination, fractional stellar radius and 
planet to star radius ratio and thus could confirm the already published results. 

Key words: planets and satellites: individual: HAT-P-18b, HAT-P-19b, HAT-P- 
27b/WASP-40b, WASP-21b 


2 M. Seeliger et al. 


1 INTRODUCTION 


Observing extra solar planets transiting their host stars has 
become an important tool for planet detection and is used to 
obtain and constrain fundamental system parameters: The 
inclination has to be close to 90°, while the planet to star 
radius ratio is constrained mainly by the transit depth. Also, 
in combination with spectroscopy, the semimajor axis and 
the absolute planet and star radii can be obtained as well. 


Several years ago, whe n the first results of the Kepler 
mission were p ublished (se e Borucki et al.l[20lll for first sci¬ 
entific results, iKoch et al.l 120101 for an instrument descrip¬ 
tion), studying the transit timing became one of the stand¬ 
ard techniques in the analysis of transit observations. Since 
space based missions are able to observe many consecut¬ 
ive transit events with a high precision, one can detect even 
small variations of the transit intervals indicating deviations 
from a strictly Keplerian motion and thus yet hidden planets 
in the observed system. Furthermore, with the discovery of 
multi-planetary systems, transit timing variations (TTVs) 
are used to find the masses of the companions without the 
need of radial velocity measurements due to the influence 
of planetary interaction on TTVs. Since many planet can¬ 
didates found in photometric surveys are too faint for radial 
velocity follow-up even with bigger telescopes, TTV analyses 
can be considered as a photometric work-around to estimate 


Although the existence of TTVs can be shown in 
already known exoplanetary systems, only a few additional 
planet candidates have been found using TTVs so far (for 
recent examples of a propo sed additional body indic ated 
by TTV analyses see e.g. iMacieiewski et ^ l2011al and 
IVan Evlen et ah] l20l4l . This is not suprising, since large 
bodies often can be found using radial velocity measure¬ 
ments or direct transit detections, while small (e.g. Earth¬ 
like) objects result in small TTV amplitudes and there¬ 
fore high precision timing measurements are needed. These 
measurements can already be acquired with medium size 
ground-based telescopes. 

Commonly the transit mid-time of each observation is 

A llotted into an observ ed minus calculated (O-C) diagram 
Ford fc HolmarJ I2OO7I ). where the difference between the 
observed transit mid-time and the mid-time obtained using 
the initial ephemeris is shown versus the observing epoch. 
In such a diagram remaining slopes indicate a wrong orbital 
period, while e.g. periodic deviations from a linear trend 
indicate perturbing forces. 

Besides the discovery of small planets, the amount of 
known massive planets on close-in orbits increased as well. 
First studies on a larger sample of planet candidates de¬ 
tected with Kepler suggest that hot giant plan ets exist in 
single plane t syste ms only (ISteffen et al.l I^12l i. However, 
ISzabo et al.l (l2013tl analysed a larger sample of Kepler hot 
Jupiters and found a few cases where TTVs can not be ex¬ 
plained by other causes (e.g. artificial sampling effects due 
to the observing cadence) but the existence of perturbers - 
additional plane ts or even exo-moon s - in the respective sys¬ 
tem. In addition ISzabo et al.l (l2013ll point toward the planet 
candidates KOI-338, KOI-94 and KOI-1241 who are all hot 
Jupiters in multi-planetary systems, as well as the HAT-P-13 
system with a hot Jupiter accompan ied by a massive p lanet 
on an eccentric outer orbit (see also lSzabo et al.ll201(]| ') and 


the WASP-12 system wi th a proposed companion candidate 
found by TTV analysis (IMacieiewski et ^12011^ '). 

The origin of those planets is yet not fully understood. 
One possible formation scenario shows that close-in giant 
planets coul d have migrated in wards after their creation 
further out JSteffen et al.l l2012h . In that case, inner and 
close outer planets would have either been thrown out of 
the system, or caught in resonance. In the latter case, even 
small perturbing masses, e.g. Earth-mass objects, can res- 
ult in TTV amplitud es in the order of seve ral minutes (see 
iFord fc HolmarLlr2007l or lSeeliger et al.l[2014h . Though Kepler 
is surveying many of those systems, it is necessary to look at 
the most promising candidates among all close-in giant plan¬ 
ets discovered so far, including stars outside the field of view 
of Kepler. In our ongoing studjQ of TTVs in exoplanetary 
systems we perform photometric follow-up observations of 
specific promising transiting planets where additional bod¬ 
ies are expected. The targets are selected by the following 
criteria: 

(i) The orbital solution of the known transiting planet 
shows non-zero eccentricity (though the circularization time- 
scale is much shorter than the system age) and/or deviant 
radial velocity (RV) data points - both possibly indicating 
a perturber. 

(ii) The brightness of the host star is U ^ 13 mag and the 
transit depth is at least 10 mmag to ensure sufficient pho¬ 
tometric and timing precision at l-2m class ground-based 
telescopes. 

(iii) The target is visible from the Northern hemisphere. 

(iv) The target has not been studied for TTV signals before. 


e xoplari ets WASP- 
2013bli . WASP- 


In the past the tra nsitini 
12b (IMacieiewski et ^ 2011a . , 

3b (IM acieiews ki et al. I2OI0I . l2013al 'l. WASP-lOb 
(|Macieiewski et al.l 201 f3. Macie jewski et al. 2014 in 
prep .), WASP-14b ( Raetd 2012ll. TrES-2 (iRaetz et al.l 
1201^ . and HAT-P-32b ( Seeliger et al.l l2014h have been 
studied by our group in detail. In most cases, except for 
WASP-12b, no TTVs could be confirmed. 

Here, we extend our investigations to search for TTVs 
in the HAT-P-18, HAT-P-19, HAT-P-27/WASP-40 and 
WASP-21 planetary systems. In Section [ 2 ] we give a short 
description of the targets analysed within this project. Sec¬ 
tion [ 3 ] explains the principles of data acquisition and reduc¬ 
tion and gives an overview of the telescopes used for obser¬ 
vation. The modeling procedures are described in Section |4l 
followed by the results in Section Finally, Section [Ogives 
a summary of our project. 


2 TARGETS 

2.1 HAT-P-18b and HAT-P-19b 

iHartman et al.l (l201ll i reported on the discovery of the exo¬ 
planets HAT-P-18b and HAT-P-19b. The two Saturn-mass 
planets orbit their early K type stars with periods of 5.51d 

and 4.Old, respectively. _ _ 

In case of HAT-P-18b IHartman et"^ (l201lh found the 


see http://ttv.astri.umk.pl/doku.php for a project overview 
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Table 1. The observing telescopes that gathered data within the TTV project for HAT-P-18b, HAT-P-19b, HAT-P-27b/WASP-40b 
and WASP-21b in order of the number of observed transit events of the Observatory. The table lists the telescopes and corresponding 
observatories, as well as the telescope diameters 0 and number of observed transit events per telescope in this project N^j.. 


# 

Observatory 

Telescope (abbreviation) 

0 (m) 

A'tr 

1 

Michael Adrian Observatory Trebur (Germany) 

TIT (Trebur 1.2m) 

1.2 

8 

2 

Graduate Institute of Astronomy Lulin (Taiwan USA) 

Tenagra II (Tenagra 0.8m) 

0.8 

5 



RCOS16 (Lulin 0.4m) 

0.4 

2 

3 

University Observatory Jena (Germany) 

90/60 Schmidt (Jena 0.6m) 

0.9/0.6 

4 



Cassegrain (Jena 0.25m) 

0.25 

3 

4 

TUBITAK National Observatory (Turkey) 

TlOO (Antalya 1.0m) 

1.0 

5 

5 

Galar Alto Astronomical Observatory (Spain) 

1.23m Telescope (CA-DLR 1.2m) 

1.23 

4 

6 

Sierra Nevada Observatory (Spain) 

Ritchey-Chretien (OSN 1.5m) 

1.5 

2 

7 

Peter van de Kamp Observatory Swarthmore (USA) 

RCOS (Swarthmore 0.6m) 

0.6 

2 

8 

National Astronomical Observatory Rozhen (Bulgaria) 

Ritchey-Chretien-Coude (Rozhen 2.0m) 

2.0 

1 



Cassegrain (Rozhen 0.6m) 

0.6 

1 

9 

Teide Observatory, Ganarian Islands (Spain) 

STELLA-1 (Stella 1.2m) 

1.2 

2 

10 

University Observatory Bochum (Gerro Armazones, Chile) 

VYSOS6 (Chile 0.15m) 

0.15 

1 

11 

Xinglong Observing Station (China) 

90/60 Schmidt (Xinglong 0.6m) 

0.9/0.6 

1 

12 

Gettysburg College Observatory (USA) 

Cassegrain (Gettysburg 0.4m) 

0.4 

1 

13 

Stara Lesna Observatory (Slovak Rep.) 

0.5m Reflector (StaraLesna 0.5m) 

0.5 

1 

14 

Istanbul University Telescope at ^anakkale (Turkey) 

0.6m Telescope (Qanakkale 0.6m) 

0.6 

1 

15 

Toruh Centre for Astronomy (Poland) 

0.6m Cassegrain Telescope (Toriin 0.6m) 

0.6 

1 


eccentricity to b e slightly non-z e ro (e = 0.084 ± 0.048). Re¬ 
cent studies of lEsposito et al.l (|2014| ) found the eccentri¬ 
city to be consistent with a non-eccentric r etrograde orbit 
by an alysing the Rossiter-McLaughlin effect. iKnutson et al.l 
ll2013l) also analysed the RV signal and found a jitter in 
the order of 17.5 ms”^ that remains unexplained. Interest¬ 
ingly, the tr ansit data listed in the exoplan et transit data¬ 
base (ETD; IPoddanv. Brat, fc Peich^ l2010l ) shows a huge 
spread in the transit depth in the order of seyeral tens of 
milli-magnitudes. Thus we included HAT-P-18b in our list 
of follow-up objects to also confirm or refute these transit 
depth variations. In addition, we performed a monitoring 
project for HAT-P-18 over a longer time span to poss ibly 
find overall brightness variations. In iGinski et al.l (l2012l l we 
used lucky imaging with Astralux at the Calar Alto 2.2m 
Telescope to search for additional low-mass stellar compan¬ 
ions in the system. With the data gathered in this previous 
study we could already exclude objects down to a mass of 
0.140±0.022 Mq at angular separations as small as 0.5 arcsec 
and objects down to 0.099±0.008 Mq outside of 2 arcsec. 

For HAT-P-19b a small eccentri city o f e = 0.067±0.042 
was determined bv iHartman et ahl (l201ll ~). They also found 
a linear trend in the RV residuals pointing towards the ex¬ 
istence of a long period perturber in the system. Within 
this project we want to address the problem of the proposed 
perturber using photometric methods, i.e. follow-up transit 
events to find planetary induced TTV signals. 


2.2 HAT-P-27b/WASP-40b 

HAT-P-27b (iBekv et al.ll201lll. mdepe ndently discovered as 
WASP- 40b by lA nderso n et al.l (l201ll l within the WASP- 
survey ijPollacco et al.ll200(ill . is a typical hot Jupiter with 
a period of 3.04d. Whi le the eccen t ricity was found to be 
e = 0.078 ± 0.047 by iBekv et ^ (l201ll ~). lAnderson et al.1 
(|201lll adopted a non-eccentric orbit. However, the latter au¬ 
thor found a huge spread in the RV data with up to 40m s“^ 
deviation from the circular single planet solution. According 


to lAnderson et al] (I2OIII ') one possible explanation, despite a 
changing activity of the K-type host star, is the exist e nce o f 
a perturber that might not be seen in the lBekv et al.l ll201ll ~) 
data due to the limited data set. However, the authors sug¬ 
gest further monitoring to clearify the nature of the system. 
One possibility is to study the companion hypothesis from 
the TTV point of view. 

Another interesting a spect of HAT-P-27b is the transit 
shape which Beky et ahl (1201 il l fi tted using a fiat bottom 
model. lAnderson et al.1 ll201ll ~) and lSada et al.l (I2OI2I I found 
the transit ra ther to have a round ish shape. From the graz¬ 
ing criterion dSmallev et al.llioill ') they concluded that the 
system is probably grazing, which would explain the unusual 
shape of t he transit. Howeve r it is not clear why this is not 
seen in the lBekv et al.1 il201ll l data, hence it is still not clear 
which shape is real. 


2.3 WASP-21b 

The planetary host star WASP-21, with its Satu r n-mas s 
planet on a 4.32d orbit discovered bv iBouchv et al.1 ll201(tl . 
is one the most metal-poor planet hosts accompanied by one 
of the least dense planets discov e red by ground-based transit 
searches to date. iBouchv et al.1 (I 2 OI 0 II found that including 
a small non-zero eccentricity to the fit does not improve the 
results. Hence, they concluded that the eccentricity is con¬ 
sistent with zero. _ _ 

However, in a later study iBarros et al.1 (|2011al l found 
the G3V star to be in the process of moving off the main 
sequence. Thus, we included further observations of WASP- 
21 b planetary transits to improve the knowledge on this sys¬ 
tem. 


3 DATA ACQUISITION AND REDUCTION 

Our observations make use of Y ETI network telescope s 
(Young Exoplanet Transit Initiativei lNeuhauser et ^l201ll l. 
a worldwide network of small to medium sized telescopes 
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Table 2. The list of all transit observations gathered within the 
TTV project sorted by object and date. Though no preselections 
for quality or completeness have been applied to this list, transits 
used for further analysis have been marked by an asterisk. The 
filter subscripts B, C and J denote the photometric systems of 
Bessel, Cousins and Johnson, respectively. The last column lists 
the number of exposures and the exposure time of each observa¬ 
tion. 


# 

Date 

Telescope 

Fiiter 

Exposures 

HAT-P-18b 

1* 

2011-04-21 

Trebur 1.2m 

Rs 

189 X 

90s 

2 

2011-05-02 

Trebur 1.2m 

Rs 

123 X 

45 s 

3* 

2011-05-24 

Trebur 1.2m 

Rs 

323 X 

60s 

4 

2011-06-04 

Rozhen 2.0m 

Rc 

1000 X 

10 s 

5* 

2012-05-05 

Rozhen 0.6m 

Ic 

219x 

90s 

6 

2012-06-07 

CA DLR 1.23m 

Bj 

250 X 

60s 

7 

2013-04-28 

Antalya 1.0m 

R 

214 X 

50 s 

8 

2014-03-30 

Toruh 0.6m 

clear 

297 X 

40 s 

HAT-P-19b 

9* 

2011-11-23 

Jena 0.6m 

Rs 

246 X 

50 s 

10 

2011-11-23 

Jena 0.25m 

Rs 

320 X 

50 s 

11 

2011-11-23 

Trebur 1.2m 

Rs 

461 X 

30 s 

12 

2011-12-05 

Jena 0.6m 

Rs 

129 X 

60s 

13 

2011-12-05 

Jena 0.25m 

Vs 

28 X 

300 s 

14* 

2011-12-09 

Jena 0.6m 

Rs 

290 X 

50 s 

15 

2011-12-09 

Jena 0.25m 

Vs 

118x 

150 s 

16 

2011-12-09 

Trebur 1.2m 

Rs 

380 X 

35 s 

17* 

2011-12-17 

CA DLR 1.23m 

Rj 

273 X 

60s 

18 

2014-08-01 

Antalya 1.0m 

R 

148 X 

60s 

19 

2014-08-05 

Antalya 1.0m 

R 

196 X 

40 s 

20 

2014-08-21 

Jena 0.6m 

Rs 

152 X 

50 s 

21* 

2014-10-04 

Jena 0.6m 

Rs 

280 X 

50 s 

HAT-P-27b 

22* 

2011-04-05 

Lulin 0.4m 

Rs 

166 X 

40 s 

23* 

2011-04-08 

Lulin 0.4m 

Rs 

250 X 

40 s 

24 

2011-05-03 

Stella 1.2m 

Ho 

180 X 

100 s 

25* 

2011-05-05 

Trebur 1.2m 

Rs 

162 X 

70s 

26 

2011-05-08 

Stella 1.2m 

Ho 

190 X 

100 s 

27 

2011-05-21 

Tenagra 0.8m 

R 

141 X 

40 s 

28 

2012-03-07 

StaraLesna 0.5m 

R 

361 X 

30 s 

29 

2012-03-29 

Tenagra 0.8m 

R 

240 X 

30 s 

30 

2012-04-01 

Tenagra 0.8m 

R 

329 X 

20 s 

31 

2012-04-04 

Tenagra 0.8m 

R 

333 X 

20 s 

32 

2012-04-25 

Xinglong 0.6m 

R 

154 X 

40 s 

33 

2012-05-16 

Trebur 1.2m 

Rs 

231 X 

70s 

34 

2012-05-25 

Chile 0.15m 

ij/Rj 

220 X 

80s 

35 

2012-06-13 

Tenagra 0.8m 

R 

223 X 

15 s 

36* 

2013-06-03 

Antalya 1.0m 

R 

156 X 

60s 

37* 

2013-06-03 

OSN 1.5m 

R 

435 X 

30 s 

38 

2013-06-06 

CA DLR 1.23m 

Rj 

172 X 

60s 

39* 

2014-06-18 

Antalya 1.0m 

R 

146 X 

50 s 

WASP-21b 

40* 

2011-08-24 

Swarthmore 0.6m 

Rs 

545 X 

45 s 

41 

2011-08-24 

Gettysburg 0.4m 

R 

230 X 

60s 

42* 

2012-08-16 

Trebur 1.2m 

Rs 

365 X 

40 s 

43 

2012-10-20 

Antalya 1.0m 

R 

242 X 

40 s 

44* 

2013-09-18 

CA DLR 1.23m 

Rj 

584 X 

30 s 

45 

2013-09-22 

Antalya 1.0m 

R 

208 X 

50 s 

46 

2013-09-22 

Ulupinar 0.6m 

Rs 

163 X 

110s 


mostly on the Northern hemisphere established to explore 
transiting planets in young open clusters. 

A summary of all participating telescopes and the num¬ 
ber of performed observations can be found in Table [T] Most 
of the observing telescopes are part of the YETI network. 
This includes telescopes at Cerro Armazones (Chile, op¬ 
erated by the University of Bochum), Gettysburg (USA), 
Jena (Germany), Lulin (Taiwan), Rozhen (Bulgaria), Sierra 
Nevada (Spain), Stara Lesna (Slovak Republic), Swarthmore 
(USA), Tenagra (USA, operated by the National Central 
University of Taiwan) and Xingl ong (China). For detail s 
about location, mirror and chip see lNeuhauser et al.l ll201lll . 

In addition to the contribution of the YETI telescopes, 
we obtained data using the following telescopes: 

• the 1.2m telescope of the German-Spanish Astronom¬ 
ical Center on Calar Alto (Spain), which is operated by Ger¬ 
man Aerospace Genter (DLR). 

• the 1.2m robotic telescope STELLA-I, situated at 
Teide Observatory on Tenerife (Spain) and operated by the 
Leibnitz-Institut fiir Astrophysik Potsdam (AIP). 

• the Trebur IMeter Telescope operated at the Michael 
Adrian Observatory Trebur (Germany) 

• the TlOO telescope of the TUBITAK National Obser¬ 
vatory (Turkey) 

• the 0.6m telescope (GIST60) at Ulupmar Observatory 
operated by Istanbul University (Turkey) 

• the 0.6m Cassegrain telescope of the Toruh Centre for 
Astronomy (Poland) 

Besides the transit observations, the Jena 
0.6m telescope with its Schmidt Teleskop Kamera 
llMugrauer fc Bertholdl I 2 OI 0 I I was used to perform a long 
term monitoring of HAT-P-18 as described in Sections 12.11 
and 15. II 

Between 2011 April and 2013 June our group observed 
45 transit events (see Table [5)) using 18 different telescopes 
(see Table [ 1 ]). 15 observations could be used for further ana¬ 
lysis, while 30 observations had to be rejected due to sev¬ 
eral reasons, e.g. no full transit event has been observed or 
bad weather cond i tions a nd hence low signal to noise. E.g. 
ISouthworth et akl (l2009al f9l showed that defocusing the tele¬ 
scope allows to reduce atmospheric and flat helding effects. 
Since a defocused image spreads the light over several CCD 
pixel, one can increase the exposure time and hence the ef¬ 
fective duty cycle of the CCD assuming a cons tant read out 
time ( as mentioned also in the conclusions of iBarros et al.1 
l2011al i . Thus we tried to defocuse the telescope and increase 
the exposure time during all our observations. Table [3] lists 
the ingress/egres s durat i ons r derived using the formula (18) 
and (19) given in IWin 3 (1201011 . With our strategy we obtain 
at least one data point within 90s. This ensures to have at 
least 10 data points during ingress/egress phase which is re¬ 
quired to ht the transit model to the data and get precise 
transit mid-times. 

All data has been reduced in a standard way W ap¬ 
plying dark/bias and flat field corrections using iraiq. The 


^ IRAF is distrubuted by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under coorporative agreement with 
the National Science Foundation. 
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Table 3. The input parameters for the JKTEBOP Sz TA P runs for all objec t s liste d in Section[2] All values have been obtained from the 
original discovery papers. LD coefficie nts are taken fromIClaret Sz BloemenI l|201lh linear interpolated in terms of ^eff’ ^^^9 and [Fe/H] 
using the EXOFAST/QUADLD code llEastman et al.ll2013[l . Free parameters are marked by an asterisk. At the bottom the duration of 
ingress and egress according to I Win 3 (l2010h has been added. 


Object 

HAT-P-15b 

HAT-P-18b 

HAT-P-19b 

HAT-P-27b 

WASP-21b 

WASP-38b 

rp -1- rs* 

0.0575(19) 

0.0575(19) 

0.0709(33) 

0.1159(65) 

0.0959(44) 

0.0829(27) 

Rp/Rs* 

0.1019(09) 

0.1365(15) 

0.1418(20) 

0.1186(31) 

0.1040(35) 

0.0844(21) 

i 

89.1(2) 

88.8(3) 

88.2(4) 

84.7(7) 

88.75(84) 

89.69(30) 

a/Rs* 

19.16(62) 

16.04(75) 

12.24(67) 

9.65(54) 

10.54(48) 

12.15(39) 

Mp/Ms 

0.0018(01) 

0.000243(26) 

0.000329(37) 

0.000663(58) 

0.000282(19) 

0.00213(11) 

e 

0.190(19) 

0.084(48) 

0.067(42) 

0.078(47) 

0 

0.0314(46) 

p(d) 

10.863502(27) 

5.5080023(06) 

4.008778(06) 

3.039586(12) 

4.322482(24) 

6.871815(45) 

R (mag) 

11.81 

12.61 

12.82 

11.98 

11.52 

9.22 

(K) 

5568(90) 

4803(80) 

4990(130) 

5300(90) 

5800(100) 

6150(80) 

logg (cgs) 

4.38(03) 

4.57(04) 

4.54(05) 

4.51(04) 

4.2(1) 

4.3(1) 

[Fe/H] (dex) 

+0.22(08) 

+0.10(08) 

+0.23(08) 

+0.29(10) 

-0.46(11) 

-0.12(07) 

V sini (km s~^) 

2.0(5) 

0.5(5) 

0.7(5) 

0.4(4) 

1.5(6) 

8.6(4) 

LD law of the star 

quadratic 

quadratic 

quadratic 

quadratic 

quadratic 

quadratic 

R band linear* 

0.4200 

0.5736 

0.5433 

0.4808 

0.3228 

0.2998 

R band non-linear* 

0.2525 

0.1474 

0.1710 

0.2128 

0.2982 

0.3095 

V band linear* 

0.5274 

0.7180 

0.6783 

0.6002 

0.4055 

0.3834 

V band non-linear* 

0.2164 

0.0697 

0.1039 

0.1643 

0.2892 

0.2992 

"^egress/ingress (mill) 

27.8 

22.8 

23.1 

37.8 

20.1 

21.9 


respective calibration images have been obtained in the same 
night and with the same focus as the scientific observations. 
This is necessary especially if the pointing of the telescope 
is not stable. When nsing calibration images obtained with 
different foci, patterns remain in the images that lead to 
distortions in the light curve. 

Besides our own observations, we also use literature 
data. This involves data from the respective d iscovery papers 
menti oned in Section as wel l as data fro m Esposito et al.l 
ll2014h in ca se of HAT-P- 1 8b. ISada et al.l (12012 1 in case of 
HAT-P-27b, iBarros e t al.| llioilal) and ICiceri et al.l (l2013l l 
for WASP-21b, and Simpson et ahl ll201oll for WASP-38b. 


4 ANALYSES 

The light curve extraction and modelling is p erformed ana¬ 
logous to the procedure described in detail in lSeelieer et al.1 
||20141 . 

4.1 Light curve extraction 

The Julian date of each image is calculated from the header 
information of the start of the exposure and the exposure 
time. To precisely determine the mid-time of the transit 
event these informations have to be stored most accurate. 
The reliability of the final light curve model thus also de¬ 
pends on a precise time synchronization of the telescope 
computer system. Observing transits with multiple tele¬ 
scopes at the same time enables to look for synchroniza¬ 
tion errors which would otherwise lead to artifacts in the 
O-C diagram (as shown for HAT-P-27b in epoch 415, see 
Section lOl) . 

We use differential aperture photometry to extract the 
light curve from the reduced images by measuring the bright¬ 
ness of all bright stars in the field with routines provided 


by IRAF. The typical aperture radius is « 1.5 times the 
mean full width half maximum of all stars in the field of 
view (FoV). The best fitting aperture is found by manually 
varying the aperture radius by a few pixels to minimize the 
photometric scatter. The final light curve is created by com¬ 
paring the star of interest against an artificial standard star 
composed of the (typically 15-30) brightest stars in the FoV 
weighted by their respective constantn ess as introduced by 
iBroeg. Fernandez, fc Neuhauserl (l2005l i. 

The final photometric errors are based on the instru¬ 
mental IRAF measurement errors. The error of the constant 
comparison stars are rescaled by their photometric scat¬ 
ter using shared scaling factors in order to achieve a mean 
Xred ~ 1 for comparison stars. The error bars of the 
transit star are rescaled afterwards using the same scaling 
facto rs (for further details on the procedure see lBroeg et al.1 

I2OO5II . 

Due to atmospheric and air-mass effects transit light 
curves show trends. They can impact the determination of 
transit parameters. To eliminate such effects we start the 
observation about 1 hour before and finish about 1 hour after 
the transit itself. Thus we can detrend the observations by 
fitting a second order polynomial to the out-of-transit data. 


4.2 Modelling with TAP and JKTEBOP 


To model the ligh t curves we used the Transit Ana¬ 
lysis Package (tap; iGazak et al.l [20l3 l. The modelling of 
the trans it light curve is done by usi ng the exofast 
routines llEastman. Gaudi. Sz Ag(3 2013ll with the light 
curve model of iMandel fc AgolT i 2002 1. For error estimation 
TAP uses Markov Chain Monte Garlo simulations (in our 
case 10 times 10® MCMC chains) together wi th wavelet- 
based likelihood functions (iGarter fc Winnll2009l i. The coef¬ 
ficients for the quadratic limb darkening (LD) law used 
by TAP are taken from the EXOFAST/QUADLD-routine of 
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lEastman et_aL (l2013lfl tha t linea rly interpolates the LD 
tables of Claret fc Bloem^ (l201lli . 

F or comparison we also use jktebop (see ISouthworthI 
I2OO8I . and references therein) which is ba s ed on the ebop 
code for eclipsing binaries dEtzell Il98ll : IPoDPer fc Etzell 
To compare the results with those obtained with 
TAP, we only use a quadratic LD law which is sufficient 
for ground-based data. For error estimation we used Monte 
Carlo simulations (10^ runs), bootstrapping (10“* data sets), 
and a residual shift method as provided by jktebop. 

As input values we use the system properties presented 
in the respective discovery papers (see Table [3] for a sum¬ 
mary). For both light curve fitting procedures we took the 
original light curve as well as a threefold binned one. Though 
the binned light curves result in a lower rms of the fit, no 
better timing result can be achieved due to a longer cadence. 

As free parameters we use the mid-transit time Tjjjjfj, 
inclination i, and planet to star radius ratio fc = rp/rs (with 
rp and rs being the planet and stellar radius scaled by the 
semimajor axis, respectively). In case of TAP the inverse frac¬ 
tional stellar radius a/Rs = 1/ts, in case of jktebop the 
sum of the fractional radii (rp -|- rs) is fitted as well. Both 
quantities are an expression of the transit duration and can 
be transformed into each other according to the following 
equation: 


a/Rs = (1 + rp/rs) / (rp -|- rs) 

The fitting procedure is applied two times. First keeping the 
limb darkening coefficients fixed at their theoretical values, 
and afterwards letting them vary. For TAP we set the fitting 
interval to ±0.2. In case of jktebop we use the option to set 
the LD coefficients fixed for the initial model, but let them 
being perturbed for the error estimation. Thus the fitted 
model does not change, but the error bars are increased. 
The eccentricity was fixed to zero for all our analyses. 

Since all data are obtained using JD^jq^Q as time base, 
we transform the fitted mid-transit times to BJDrpjqg after¬ 
wards using the online converteiQ provided by Jason East¬ 
man (for a detailed description of the bary centric dynamical 
time see lEastman. Siverd. fc GaudiluOlOl 'l. 

Finally, we d erive the photometric noise rate (pnr, 
iFulton et al.ll201ll ~) as a quality marker for all light curves, 
which is defined as the ratio between the root mean square of 
the model fit and the number of data points per minute. For 
further analysis we took data with pnr < 4.5 into account. 


5 RESULTS 

For every light curve we get six different models, four from 
TAP (for the binned and unbinned data with the LD coeffi¬ 
cients fixed and free) and two from jktebop (for the binned 
and unbinned data, LD coefficients set free for error estim¬ 
ation only). To get one final result we averaged those six 
results. As for the errors, we got four different estimations 
from TAP and 12 from jktebop. As final error value, we took 
the maximum of either the largest of the error estimates, or 


® the limb darkening calculator is available online at 
http: //astroutils.astronomy.ohio-state.edu/exofast/limbdark.shtml 
^ http://astroutils.astronomy.ohio-state.edu/time/utc2bjd.html 


HAT-P-18 brightness monitoring 



Figure 1. Relative R band brightness of the star HAT-P-18 over 
a time span of 12 months. The dotted line represents the rms of 
a constant fit. 


the spread of the model fit results to use a conservative er¬ 
ror estimate. It has to be noted, though, that the spread 
between the different models has always been below size of 
the error bars. 

Binning the light curve threefold using an error 
weighted mean in principle still leaves enough data points 
during ingress and egress to be able to fit the transit 
model to the data while reducing the error bar of an in¬ 
dividual measurement. However, comparing the results of 
the threefold binned and unbinned data we do not see sig¬ 
nificant differences, neither in the fitted values, nor in the 
error bars. 

The same counts for the differences between the TAP 
models obtained with fixed LD values and those obtained 
with the LD coefficients set free to fit. For a detailed discus¬ 
sion of the influence of the LD model on transit light curves 
see e.g. iRaetz et al.l ll2ni4h . 


5.1 HAT-P-18b 

For HAT-P-18b we obtained five transit light curves (see 
Fig [ 2 ) 1 . All light curves show some features which could 
be caused by stellar activity, e.g. spots. However, the qual¬ 
ity of the data does not allow to draw further conclusion. 
Moreover, there is only a small number of suitable compar¬ 
ison stars available in the respective FoV of each observation, 
hence this could also be an artificial effect. 

Except for one - but not significant - outlier the differ¬ 
ences in the O-C diagram (see Fig. [2]) can be explained by 
redetermining the published period by (0.53 ± 0.36) s. Hence 
we find a slightly l arger period compare d to the originally 
proposed period of lHartman et aP ll201ll '>. 

Regarding the transit depth, and thus the planet to star 
radius ratio, we do not see any changes. However, having a 
look at the data provided in the exopl anet transit database 
fETD: [Poddanv. Brat, fc Peichall2010ll one can see that the 
values for the transit depth reported there vary by several 
tens of mmag. Such transit depth variations can be caused 
by close variable stellar companions placed within the aper¬ 
ture due to the pixel scale of our detectors and the defocus¬ 
sing of the telescopes. Thus we took a more detailed look at 
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epoch 174- 2011-04-21 - HAT-P-18b @ Trebur 1.2m 



epoch 180 - 2011-05-24 - HAT-P-18b @ Trebur 1.2m 



time [JD-2455706) 


epoch 243-2012-05-05 - HAT-P-18b @ Rozhen 0.6m 
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Figure 2. top: The threefold binned transit light curves of the three complete transit observations of HAT-P-18b. The upper panels 
show the light curve, the lower panels show the residuals. The rms of the fit of the threefold binned light curves (dotted lines) are shown 
as well, bottom: The present result for the HAT-P-18b observing campaign, including the O-C diagram, as well as the results for the 
rever se fractional ste l lar ra dius a/Rs, the inclination i, and the planet to star radius ratio k. The open circle denotes literature data 
from [Hartman et al.l ll201lh , the open triagles denotes data from Esposito et al.l ||2014| V Filled triangles denote our data (from Trebur 
and Rozhen. The dotted line shows the Icr error bar of the constant fit. 


the long term variability of the parent star. Over a timespan 
of 12 months we obtained 3 images in four bands (B, V, R, I) 
in each clear night using the Jena 0.6m telescope. As shown 
in Fig. m the mean variation of the R band brightness is 
~ 3.8 mmag taking the individual error bar of the meas¬ 
urements into account. However, this variation is too small 
to be responsible for the seen transit depth variations. The 
monitoring of the remaining bands is not shown but leads to 
a similar result. Having a closer look at the light curves listed 
for HAT-P-18b in ETD, one can see that a large number is 
of lower quality. This is especially true for those light curves 
responsible for the spread in the tabulated transit depth. 
Taking only the higher quality data into account the spread 
is much smaller. Depending on the quality cut the variation 
can even reach the order of the error bars of the measure¬ 


ments. Thus we believe that the transit depth variation is 
negligible. 

Despite a spread in the data, which can be explained by 
the quality of the light curves, we do not see any significant 
difference for k, i and a/f?s between the respective observa¬ 
tions. A summary of all obtained parameters can be found 
in Table [4] at the end of the paper, as well as a comparison 
with literature values. 

5.2 HAT-P-19b 

For HAT-P-19b we got two light curves using the Jena 
0.6m and one light curve from the CAHA 1.2m telescope 
(Fig. [ 3 ]). In all three cases we obtained high precision 
data. The light curves show no artifacts that could be 
ascribed to e.g. spots on the stellar surface. Plotting the 
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mid-transit times into the O-C diagram, we can redeter¬ 
mine the period by (0.53 ± 0.06) s. As for the inclination 
and the reverse fr actional stellar r a dius w e can conhrm the 
values reported in iHartman et aP (1201 ll i. The radius ratio 
of fc = 0.13 78 ± 0.0014, however, seems to be smaller than 
assumed bv IHartman et ^ (I 2 OIII I (fe = 0.1418 ± 0.0020). 


5.3 HAT-P-27b 

HAT-P-27b planetary transits were observed six times. An 
advantage of a network such as YETI lies within the possib¬ 
ility of simultaneous observations using different telescopes. 
This enables us to independently check whether the data is 
reliable. For HAT-P-27b simultaneous observations could be 
achieved at epoch 415 using two different telescopes (Anta¬ 
lya 1.0m and OSN 1.5m). 

As one can see in Fig. |T] the best-htting transit shapes 
differ. Since some data has been acquired using small tele¬ 
scopes under unfavorable conditions, this might be an ar¬ 
tificial effect. However, these sha pe variations can a lso be 
seen in the literature data. While iBekv et al.l 1I2OIIII fitted 
a flat bottom to the transit they observed, the transit of 
lAnderson etaljJ gOllIl shows a rather roun dish bottom and 
so doe s the Sadu^_et_^ (|201^ transit. Both [Anderson et al.l 
(120111) and ISada et al.l (I2OI2I) claim, a roundish bottom is 
in good agreement with a grazing transit. Our most precise 
transit light curve of HAT-P-27b at epoch 415 the OSN 1.5m 
telescope also shows no flat bottom. Thus, we would agree 
with the previous authors that the planet is grazing. The 
epoch 540 observations also shows a V-shape. However, due 
to a connection problem parts of the ingress data is missing, 
thus precise fits of the transit parameters are not possible. 

In addition to o ur own data we a dded data from 
ISada et al.l (I2OI2I ) and [Brown et al.l (I2OI2I) . The latter one 
only lists system parameters without giving an epoch of ob¬ 
servation, thus we artihcially put them to epoch 200. Since 
we do not see any systematic trend in the remaining data, 
this does not effect the conclusion on the system parameters 
but improves the precision of the constant fit. 

The system parameters i, a/Rs and k can be determ¬ 
ined more precisely than before taking the errors of the in¬ 
dividual measurements into account. All three parameters 
are in good agreement with the results of previous authors. 
Furthermore, we do not see any significant variation. The 
larger fc-value of the epoch 540 observations are due to the 
quality of the corresponding light curve. 

Looking at the mid-transit time we see that a period 
change of (—0.51 ± 0.12) s explains the data quite well. The 
mid-transit time of one of the epoch 415 observations was 
found to be « 4.5 min ahead of time, while the other one is 
as predicted. This way we could identify a synchronization 
error during one of the observations. This example shows 
the importance of simultaneous transit observations. Unfor¬ 
tunately this was the only sucessful observation of that kind 
within this project ( for a larger set of d ouble and threefold 
observations see e.g. [Seelieer et aL[[2014l ). 


5.4 WASP-21b 

Four transit light curves of W ASP-21b are availab le, in¬ 
cluding one light curve from [Barros et al.[ ([2011al) (see 


FigEj). In addition, the resu lts of the analysis of two transit 
events of Ciceri_etal. 1(H) and one transit observation of 
ISouthworthl ( 20121) are also taken into account. Concern¬ 
ing the O-C diagram, we found that a period change of 
(2.63 ± 0.17) s removes the linear trend which is present in 
the data fitted with the initial ephemeris. As in the previ¬ 
ous analyses no trend or sinusoidal variation in the system 
parameters can be seen. 

However, regarding inclination and reverse fractional 
stellar radius we do see a signihcant differ ence between 
our re sults and the initial values published bv [Bouchv et al.[ 
(I2OIOI) . T his was also found by other authors before. As dis¬ 
cussed in [Barros et al. |[2011a[) this result is a consequence 
of the assumption of Bouchv et al.l (l201(j) that the planet 
host star is a main sequence star, while [Barros et ^ (l2011ah 
found that the star starts evolving off the main sequence and 
thus its radius increases. This in turn leads to corrections of 
the stellar and hence planetary properties. 


6 SUMMARY 

We presented the results of the transit observations of 
the extra solar planets HAT-P-18b, HAT-P-19b, HAT-P- 
27b/WASP-40b and WASP-21b which are part of our on¬ 
going project on gound-based follow-up observations of exo¬ 
planetary transits using small to medium sized telescopes 
with the help of YETI network telescopes. During the past 
three years we followed these well chosen objects to refine 
their orbital parameters as well as to find transit timing 
variations indicating yet unknown planetary companions. 
Table H contains an overview of the redetermined proper¬ 
ties, as well as the available literature values, while Table [5] 
lists the results of the individual light curve fits. 

In all cases we could redetermine the orbital paramet¬ 
ers. Especially the period could be determined more precise 
than before. So far, we can not rule out the existence of TTV 
signals for the planets investigated within this study due to 
the limited number of available high quality data. Also the 
parameters a/Rs, rp/rs and inclination have been obtained 
and compared to the available literature data. Despite some 
corrections to the literature data, we found no signihcant 
variations within these parameters. To distinguish between 
a real astrophysical source of the remaining scatter and ran¬ 
dom noise as a result of the quality of our data more high 
precision transit observations would be needed. 

HAT-P-18b was also part of an out-of-transit monitor¬ 
ing for a spread in the transit depth was reported in the 
literature that could be due to a signihcant variability of 
the transit host star. Regarding our transit data we can not 
conhrm the spread in transit depth. Looking at the quality 
of the literature data showing the transit depth variation, it 
is very likely that this spread is of artihcial nature. Thus it 
is not suprising that we did not hnd stellar variability larger 
than « 3.8mmag. However we do see some structures in the 
light curves that could be caused by spot activity on the 
stellar surface. 
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Figure 3. top: The transit light curves obtained for HAT-P-19b. bottom: The present result fo r the HAT-P-19b observing campaign. 
All explanations are equal to Fig. [2] The open circle denotes literature data from I Hartman et al.l (120111) . filled triangles denote our data 
(from Jena and Calar Alto). 
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Figure 5. top: The transit light curves obtained for WASP-21b. bottom: The present result for the WASP-21b o bserving campaign. All 
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triangles denote our data (from Swarthmore, Trebur and Calar Alto). 
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Table 4. A comparison between the results obtained in Our analysis and the literature data. All epochs To are converted to BJDj'db- 


To{d) 


p(d) 


a/Rs 


k = Rp / Rs 


in 


HAT-P-18b 


our analysis 

Hartman et al. ('2011') 
EsDosito et al. ('2014') 

2 454 715.022 54 ±0.000 39 

2 454 715.022 51 ± 0.000 20 

2 455 706.7 ±0.7 

5.508 0291 ±0.000 004 2 
5.508 023 ±0.000 006 
5.507 978 ±0.000 043 

17.09 ±0.71 
16.04 ±0.75 
16.76 ±0.82 

0.136 2 ±0.0011 
0.136 5 ±0.001 5 
0.136 ±0.011 

88.79 ±0.21 
88.3 ±0.3 
88.79 ±0.25 

HAT-P-19b 

our analysis 

Hartman et al. ('2011') 

2 455 091.535 00 ±0.000 15 

2 455 091.534 94 ±0.000 34 

4.008 784 2 ±0.000 000 7 
4.008 778 ±0.000 006 

12.36 ±0.09 
12.24 ±0.67 

0.1378±0.0014 
0.141 8 ±0.002 0 

88.51 ±0.22 
88.2 ±0.4 

HAT-P-27b 

our analysis 

Bekv et al. ('2011) 
Anderson et al. ('2011) 
Sada et al. (2012) 
Brown et al. (2012) 

2 455 186.019 91 ±0.000 44 

2 455 186.019 55 ±0.000 54 

2 455 368.394 76 ±0.000 18 

2 455 186.198 22 ± 0.000 32 

3.039 580 3 ±0.000 0015 

3.039 486 ±0.000 012 

3.039 5721 ±0.000 007 8 

3.039 582 4 ±0.000 003 5 

3.039 577 ±0.000 006 

10.01 ±0.13 

y.oo _Q 4Q 

9.88 ±0.39 

9-11 

9.80 

0.119 2 ±0.001 5 

0.118 6 ±0.0031 

0.125 0 ±0.001 5 

u.ro‘±‘± -0.038 9 

0 120 +0009 

U.IZU -0,007 

85.08 ±0.07 
84.7 t°;I 
84.981°;™ 
84.23 ±0.88 

85.0 ±0.2 


WASP-21b 


our analysis 
^ouch^^t_^ (20W) 
^arros_et_^^ (201^) 
(20^) 

^outhwor^ (20^) 


2 454 743.042 17 ± 0.000 65 
2 454 743.042 6 ±0.002 2 
2 455 084.520 48 ±0.000 20 
2 454 743.040 54 ±0.000 71 
2 455 084.520 40 ±0.000 16 


4.322 512 6 ±0.000 002 2 

A “^99 zlH9 +0.000 024 
_o.000 019 

4.322 506 0 ±0.000 0031 
4.322 518 6 ±0.000 003 0 
4.322 506 0 ±0.000 0031 


9.62±0.17 
6.05l“:°i 
9.681°;?° 
9.46 ±0.27 
9.35 ±0.34 


0.1030±0.0008 87.12±0.24 


0 104 0 

U.IU+U _0,0018 

0 107 1 ® 

U.iU/ 1 _0,0008 

0.1055 ±0.0023 
0.1095 ±0.0013 


88 . 751 °;™ 
87.34 ±0.29 
86.97 ±0.33 
86.77 ±0.45 


Table 5. The results of the induvidual fits of the observed complete transit event. The rms of the fit and the resultant pnr are given in 
the last column. The table also shows the result for the transits with pnr > 4.5 that are not used for redetermining the system properties. 


date 

epoch 

telescope 

Tjjjid - 2 450 000 d 

a/Rs 

k = Rp / Rs 

in 

rms 1 pnr (mmag) 

HAT-P-18b 

2011/04/21 

174 

Trebur 1.2m 

5 673.419 67 ±0.001 24 

16.4 ±1.4 

0.139 9 ±0.007 2 

88.52 ± 0.84 

3.0/3.3 

2011/05/24 

180 

Trebur 1.2m 

5 706.469 93 ±0.000 80 

18.28 ±0.83 

0.134 3 ±0.003 9 

89.52 ±0.58 

3.7/4.0 

2012/05/05 

243 

Rozhen 0.6m 

6 053.472 76 ±0.000 84 

16.04 ± 1.36 

0.137 3 ±0.004 7 

88.55 ±0.79 

3.5/4.4 

2012/06/07 

249 

CA-DLR 1.2m 

6 086.518 56 ±0.001 25 

- 

- 

- 

4.1/4.5 

2013/04/28 

308 

Antalya 1.0m 

6 411.496 38 ±0.000 84 

15.22 ± 1.52 

0.1464±0.0068 

87.89 ± 0.75 

3.9/5.0 

HAT-P-19b 

2011/11/23 

199 

Jena 0.6m 

5 899.283 45 ±0.000 49 

12.56 ±0.34 

0.136 9 ±0.002 6 

88.20 ± 0.64 

2 .1/2.2 

2011/12/09 

203 

Jena 0.6m 

5 905.318 10 ±0.000 44 

12.29 ±0.35 

0.136 9 ±0.002 3 

89.05 ±0.67 

2.3/2.4 

2011/12/17 

205 

CA-DLR 1.2m 

5 913.335 71 ±0.000 34 

11.96 ±0.53 

0.136 8 ±0.002 7 

88.38 ± 0.80 

1.2/1.3 

2014/10/04 

460 

Jena 0.6m 

6 935.575 59 ±0.000 55 

12.43 ±0.36 

0.134 0 ±0.002 6 

89.25 ±0.67 

2 .8/3.0 

HAT-P-27b 

2011/04/05 

155 

Lulin 0.4m 

5 657.153 33 ±0.001 07 

10.72 ± 1.67 

0.123 3 ±0.0081 

85.53 ±0.93 

3.4/3.8 

2011/04/08 

156 

Lulin 0.4m 

5660.19481 ±0.001 16 

9.43 ± 1.01 

0.1228±0.0149 

84.69 ±0.81 

3.4/3.2 

2011/05/05 

165 

Trebur 1.2m 

5 687.551 22 ±0.000 51 

9.83 ±0.56 

0.115 3 ±0.002 9 

85.07 ±0.40 

1 .6/1.8 

2012/04/01 

274 

Tenagra 0.8m 

6 018.864 57 ±0.002 32 

9.65 ± 1.63 

0.1199±0.0126 

84.13 ±1.63 

5.7/5.8 

2012/04/25 

282 

Xinglong 0.6m 

6 043.180 95 ±0.001 35 

9.89 ± 1.67 

0.118 6 ±0.006 7 

84.83 ± 1.24 

4.3/5.1 

2013/06/03 

415 

Antalya 1.0m 

6 447.442 68 ±0.001 66 

10.64 ± 1.30 

0.1184 ±0.0081 

85.51 ±0.94 

2.6/3.5 

2013/06/03 

415 

OSN 1.5m 

6 447.445 71 ±0.000 30 

10.18 ±0.29 

0.122 4 ±0.003 7 

85.23 ±0.21 

1.2/0.9 

2013/06/18 

540 

Antalya 1.0m 

6 827.395 45 ±0.002 20 

10.77 ± 1.01 

0.146 2 ±0.0141 

85.26 ±0.55 

3.1/4.0 

WASP-21b 

2011/08/24 

244 

Swarthmore 0.6m 

5 797.734 00 ±0.001 12 

9.94 ±0.93 

0.101 4 ±0.003 2 

87.74 ± 1.41 

3.3/3.1 

2012/08/16 

327 

Trebur 1.2m 

6156.502 60 ±0.001 15 

9.97 ±0.92 

0.101 7 ±0.003 2 

87.78 ± 1.46 

2.9/2.6 

2013/09/18 

420 

CA-DLR 1.2m 

6 558.496 48 ±0.000 73 

9.38 ±0.69 

0.1064 ±0.002 7 

86.91 ±0.96 

1.6/1.4 
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